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What is claimed is: 

1. A method for estimating a seismic velocity field from seismic data 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations, the 
method comprising: 

a) defining a control plane having an edge intersecting a plurality of 
the CMP locations; 

b) producing an initial velocity field for said control plane, said initial 
velocity field comprising a plurality of time-velocity values for 
each of the CMP locations; and 



c) producing an optimized velocity field for said control plane by 
adjusting said time-velocity values for each of the CMP locations 
in response to trends, relative to offset distance, in time values, 
associated with common seismic events, until said optimized 
velocity field satisfies a condition. 

2. The method of claim 1 wherein producing said optimized velocity field 
comprises producing a set of gathers in response to said initial velocity 
field, according to a type of velocity field desired to be produced. 

3. The method of claim 2 wherein producing said set of gathers 
comprises producing sets of common offset migration image gathers in 
response to the seismic data and said initial velocity field when an 
imaging velocity field is desired to be produced. 
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The method of claim 2 wherein producing said set of gathers 
comprises performing a normal moveout operation on the seismic data 
when said imaging velocity field is not desired to be produced. 

The method of claim 2 further comprising producing an adjusted 
velocity field in response to said set of gathers and said initial velocity 
field. 

The method of claim 5 wherein producing said adjusted velocity field 
comprises, for said each CMP location and for each of said common 
seismic events, finding a slope of a curve approximating a trend, 
relative to offset, in said time values, associated with said each 
common seismic event. 

The method of claim 6 wherein finding said slope of said curve 
comprises producing slant stack values for a plurality of time(tau)-slope 
pairs. 

The method of claim 7 further comprising, for said each CMP location, 
producing for each time value of said time(tau)-slope pairs a best slope 
value in response to said time(tau)-slope pairs and said slant stack 
values. 

The method of claim 8 further comprising producing a temporary time- 
velocity representation in response to said initial velocity field and a 
plurality of said best slope values. 

The method of claim 9 wherein producing said temporary time-velocity 
representation comprises adjusting each velocity value at a given time t 
of a corresponding time-velocity representation of said initial velocity 
field according to a relation: 
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Vadj& = 




is an adjusted velocity value associated with said given 
time t\ 

is a Tau value equal to said given time /; 
is said best slope value at r = t ; and 

is said each velocity value at said given time t of said 

corresponding time-velocity representation. 

The method of claim 10 further comprising conditioning said plurality of 
best slope values to produce conditioned best slope values and 
producing an adjusted time-velocity representation in response to said 
temporary time-velocity representation and said conditioned best slope 
values according to said relation. 

The method of claim 5 further comprising determining whether said 
adjusted velocity field meets optimization criteria. 

The method of claim 12 wherein determining whether said adjusted 
velocity field meets optimization criteria comprises determining whether 
velocity values of said adjusted velocity field are within a range of 
corresponding velocity values of a temporary initial velocity field. 



where 

VadjV 

T 

P T 

V rep 
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14. The method of claim 12 further comprising producing a new adjusted 
velocity field using said adjusted velocity field as a new temporary 
initial velocity field when said optimization criteria are not met. 

5 15. The method of claim 12 further comprising identifying said adjusted 

velocity field as said optimized velocity field for said control plane when 
said optimization criteria are met. 

16. The method of claim 1 wherein producing said initial velocity field 
10 comprises producing a starting velocity field estimate. 

17. The method of claim 16 wherein producing said starting velocity field 
estimate comprises defining a range of velocity values, defining one or 
more ranges of time values and associating each velocity value of said 

15 range of velocity values with one or more corresponding ranges of time 
values. 

18. The method of claim 17 further comprising, for each of a plurality of 
CMP gathers, defining a window in which a selected CMP location 

20 associated with said each CMP gather is centered. 

19. The method of claim 18 further comprising producing a semblance 
panel associated with said selected CMP location, in response to said 
range of velocity values and selected CMP gathers associated with 

25 CMP locations within said window. 



20. The method of claim 19 further comprising producing a time-velocity 
profile associated with said selected CMP location in response to said 
semblance panel and a set of velocity stacking weights. 

30 

21. The method of claim 20 further comprising producing smoothing 
weights for each velocity value of said range of velocity values and 
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each time value of said one or more correspondence ranges of time 
values in response to respective products of total semblance and 
reciprocal velocity gradient associated with said each time value and 
said each velocity value. 

5 

22. The method of claim 21 further comprising producing a smooth time- 
velocity profile associated with said selected CMP location in response 
to said time-velocity profile and said smoothing weights. 

10 23. The method of claim 22 further comprising producing laterally edited 

and laterally smoothed time-velocity profiles in response to a plurality 
of said smooth time-velocity profiles. 

24. The method of claim 16 further comprising producing a migrated 
1 5 starting velocity field. 

25. The method of claim 24 wherein producing said migrated starting 
velocity field comprises, for said each CMP location producing a 
normal moveout gather in response to said starting velocity field 

20 estimate and seismic data associated with said each CMP location. 

26. The method of claim 25 further comprising producing a CMP stack for 
said each CMP location in response to a corresponding said normal 
moveout gather. 

25 

27. The method of claim 26 further comprising producing a simulated CMP 
gather for said each CMP location in response to a corresponding said 
CMP stack and a corresponding time-velocity profile associated with 
said each CMP location. 
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28. The method of claim 27 further comprising producing pre-stack imaged 
simulated gathers in response to a plurality of respective said 
simulated CMP gathers and respective said CMP locations. 



5 29. The method of claim 28 further comprising producing normal moveout 

pre-stack imaged simulated gathers in response to said pre-stack 
imaged simulated gathers, said CMP locations and said starting 
velocity field estimate. 

10 30. The method of claim 29 wherein producing said migrated velocity 

comprises producing said migrated velocity field in response to said 
normal moveout pre-stack imaged simulated gathers and said starting 
velocity field estimate. 



15 31. The method of claim 30 wherein producing said migrated velocity field 

comprises, for said each CMP location and for each of said common 
seismic events, finding a slope of a curve approximating a trend, 
relative to offset, in said time values, associated with said each 
common seismic event. 

20 

32. The method of claim 31 wherein finding said slope of said curve 
comprises producing slant stack values for a plurality of time(tau)-slope 
pairs. 

25 33. The method of claim 32 further comprising, for said each CMP location, 

producing for each time value of said time(tau)-slope pairs a best slope 
value in response to said time(tau)-slope pairs and said slant stack 
values. 



30 



34. 



The method of claim 33 further comprising producing a first temporary 
time-velocity representation in response to said starting velocity field 
estimate and a plurality of said best slope values. 
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35. The method of claim 34 wherein producing said first temporary time- 
velocity representation comprises adjusting each velocity value at a 
given time t of a corresponding time-velocity representation of said 
starting velocity field estimate according to the relation: 



where 



10 Vadj^ is an ac U usted velocity value associated with said given 

time /; 

r is a Tau value equal to said given time 

p T is said best slope value at r = / ; and 

V rep is sa ' d each velocity value at said given time t of said 

1 5 corresponding time-velocity representation. 

36. The method of claim 35 further comprising conditioning said plurality of 
best slope values to produce conditioned best slope values and 
producing an adjusted time-velocity representation in response to said 
20 temporary time-velocity representation and said conditioned best slope 

values according to said relation. 



25 



37. 



The method of claim 24 further comprising producing pre-stack imaged 
gathers by performing a 2-dimensional pre-stack imaging operation on 
seismic data associated with said control plane. 
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38. The method of claim 37 further comprising producing normal moveout 
gathers in response to said migrated starting velocity field by 
performing respective normal moveout operations on said pre-stack 
imaged gathers. 

5 

39. The method of claim 38 further comprising adjusting said migrated 
starting velocity field in response to said normal moveout gathers and 
said migrated starting velocity field. 

10 40. The method of claim 39 wherein adjusting said migrated starting 

velocity field comprises, for said each CMP location and for each of 
said common seismic events, finding a slope of a curve approximating 
a trend, relative to offset, in said time values, associated with said each 
common seismic event. 

15 

41. The method -of -claim 40 wherein finding said slope of said curve- 
comprises producing slant stack values for a plurality of time(tau)-slope 
pairs. 

20 42. The method of claim 41 further comprising, for said each CMP location, 

producing for each time value of said time(tau)-slope pairs a best slope 
value in response to said time(tau)-slope pairs and said slant stack 
values. 

25 43. The method of claim 42 further comprising producing a second 

temporary time-velocity representation in response to said migrated 
starting velocity field and a plurality of said best slope values. 



44. The method of claim 43 wherein producing said second temporary 
30 time-velocity representation comprises adjusting each velocity value at 

a given time t of a corresponding time-velocity representation of said 
migrated starting velocity field according to the relation: 
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Vadj(0 = 



*Pr + — 

V rep 



where 



10 



VadjiO is an adjusted velocity value associated with said given 



T 

P T 

V rep 



time t\ 

is a Tau value equal to said given time t ; 
is said best slope value at r = t ; and 

is said each velocity value at said given time t of said 

corresponding time-velocity representation. 



15 



45. The method of claim 44 further comprising conditioning said plurality of 
best slope values to produce conditioned best slope values and 
adjusting said migrated starting velocity field in response to said 
temporary time-velocity representation and said conditioned best slope 
values according to said relation. 



20 



46. The method of claim 1 further comprising determining whether the 
seismic data is 2-dimensional data or 3-dimensional data, determining 
whether an imaging velocity is desired to be produced, and performing 
a directional migration on the seismic data when the seismic data is 3- 
dimensional data and an imaging velocity is desired to be produced. 



25 47. The method of claim 46 further comprising determining whether or not 

a dip moveout velocity or said imaging velocity is desired to be 
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produced, and if so, performing a dip moveout operation on the seismic 
data. 

The method of claim 47 wherein performing said dip moveout operation 
comprises performing a 2-dimensional dip moveout operation when the 
seismic data is said 2-dimensional seismic data and performing a 3- 
dimensional dip moveout operation when the seismic data is 3- 
dimensional data. 

A method for estimating a seismic velocity field from seismic data 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations 
aligned along edges of respective planes in a 3-dimensional region in 
the earth, the method comprising: 

a) designating one of said respective planes as a control plane; 

b) estimating a velocity field for said control plane, according to the 
method of claim 1; and 

c) producing an optimized velocity field for one of said respective 
planes nearby said control plane by adjusting time-velocity 
values for each of the CMP locations associated with said 
control plane, in response to trends, relative to offset, in time 
values, associated with common seismic events, in seismic 
data associated with said one of said planes nearby said control 
plane until said optimized velocity field satisfies a condition. 

A method for estimating a seismic velocity field from seismic data 
comprising time-amplitude representations associated with source- 
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receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations 
aligned along edges of respective planes in a 3-dimensional region in 
the earth, the method comprising: 

a) designating one of said respective planes as a control plane; 

b) estimating a velocity field for said control plane, according to the 
method of claim 1 ; and 

c) producing an optimized velocity field for each of the respective 
planes in the 3-dimensional region by adjusting time-velocity 
values for each of the CMP locations associated with said 
control plane, in response to trends, relative to offset, in time 

— - values, associated with common seismic events, relative to 
offset, in the seismic data associated with said each of said 
respective planes until said optimized velocity field satisfies a 
condition. 

A method for estimating a seismic velocity field from seismic data 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations 
aligned along edges of respective planes in a 3-dimensional region in 
the earth, the method comprising: 

a) designating one of the respective planes as a temporary control 
plane; 
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b) estimating a velocity field for said temporary control plane, 
according to the method of claim 1; and 

c) producing an optimized velocity field for successive ones of the 
respective planes in the 3-dimensional region by, for a given 
plane in the 3-dimensional region, adjusting time-velocity values 
of a predecessor one of the respective planes, starting with said 
temporary control plane, in response to trends, relative to offset, 
in time values, associated with common seismic events and a 
common CMP, in seismic data associated with said given plane 
until said optimized velocity field satisfies a condition. 

A computer readable medium encoded with codes for directing a 
processor circuit to carry out the method of claim 1. 

A computer readable signal encoded with codes for directing a 
processor circuit to carry out the method of claim 1. 

An apparatus for estimating a seismic velocity field from seismic data 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations, the 
apparatus comprising: 

a) means for defining a control plane having an edge intersecting a 
plurality of the CMP locations; 

b) means for producing an initial velocity field for said control 
plane, said initial velocity field comprising a plurality of time- 
velocity values for each of the CMP locations; and 
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c) means for producing an optimized velocity field for said control 
plane by adjusting said time-velocity values for each of the CMP 
locations in response to trends, relative to offset distance, in 
time values, associated with common seismic events, until said 
optimized velocity field satisfies a condition. 

An apparatus for estimating a seismic velocity field from seismic data 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations, the 
apparatus comprising: 

a) a component configured to define a control plane having an 
edge intersecting a plurality of the CMP locations; 



b) a component configured to produce an initial velocity field for 
said control plane, said initial velocity field comprising a plurality 
of time-velocity values for each of the CMP locations; and 

c) a component configured to produce an optimized velocity field 
for said control plane by adjusting said time-velocity values for 
each of the CMP locations in response to trends, relative to 
offset distance, in time values, associated with common seismic 
events, until said optimized velocity field satisfies a condition. 

The apparatus of claim 55 wherein said component configured to 
produce said optimized velocity field is configured to produce a set of 
gathers in response to said initial velocity field, according to a type of 
velocity field desired to be produced. 
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57. The apparatus of claim 56 wherein said component configured to 
produce said optimized velocity field is configured to produce said set 
of gathers by producing sets of common offset migration image gathers 
in response to the seismic data and said initial velocity field when an 

5 imaging velocity field is desired to be produced. 

58. The apparatus of claim 56 wherein said component configured to 
produce said set of gathers by performing a normal moveout operation 
on the seismic data when said imaging velocity field is not desired to 

10 be produced. 

59. The apparatus of claim 56 further comprising a component configured 
to produce an adjusted velocity field in response to said set of gathers 
and said initial velocity field. 

15 

60. The apparatus of claim 59 wherein said component configured to 
produce said adjusted velocity field comprises, a processor circuit 
configured to: for said each CMP location and for each of said common 
seismic events, find a slope of a curve approximating a trend, relative 

20 to offset, in said time values, associated with said each common 

seismic event. 

61 . The apparatus of claim 60 wherein said processor circuit is configured 
to find said slope of said curve by producing slant stack values for a 

25 plurality of time(tau)-slope pairs. 

62. The apparatus of claim 61 wherein said processor circuit is configured 
to, for said each CMP location, produce for each time value of said 
time(tau)-slope pairs a best slope value in response to said time(tau)- 

30 slope pairs and said slant stack values. 
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The apparatus of claim 62 wherein said processor circuit is configured 
to produce a temporary time-velocity representation in response to said 
initial velocity field and a plurality of said best slope values. 



The apparatus of claim 63 wherein said processor circuit is configured 
to produce said temporary time-velocity representation by adjusting 
each velocity value at a given time t of a corresponding time-velocity 
representation of said initial velocity field according to a relation: 



1 



V 



V rep 



where 



VadjiO 

T 

P T 

Vrep 



is an adjusted velocity value associated with said given 
time /; 

is a Tau value equal to said given time t \ 
is said best slope value at r = t ; and 

is said each velocity value at said given time t of said 

corresponding time-velocity representation. 



The apparatus of claim 64 wherein said processor circuit is configured 
to condition said plurality of best slope values to produce conditioned 
best slope values and to produce an adjusted time-velocity 
representation in response to said temporary time-velocity 
representation and said conditioned best slope values according to 
said relation. 
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66. The apparatus of claim 59 further comprising a processor circuit 
configured to determine whether said adjusted velocity field meets 
optimization criteria. 

5 67. The apparatus of claim 66 wherein said processor circuit is configured 

to determine whether said adjusted velocity field meets optimization 
criteria by determining whether velocity values of said adjusted velocity 
field are within a range of corresponding velocity values of a temporary 
initial velocity field. 

10 

68. The apparatus of claim 66 wherein said processor circuit is configured 
to produce a new adjusted velocity field using said adjusted velocity 
field as a new temporary initial velocity field when said optimization 
criteria are not met. 

15 

-69. The apparatus of claim 66 wherein said processor circuit is configured 
to identify said adjusted velocity field as said optimized velocity field for 
said control plane when said optimization criteria are met. 

20 70. The apparatus of claim 55 wherein said component configured to 

produce said initial velocity field comprises a processor circuit 
configured to produce a starting velocity field estimate. 

71. The apparatus of claim 70 wherein said processor circuit is configured 
25 to produce said starting velocity field estimate by defining a range of 

velocity values, defining one or more ranges of time values and 
associating each velocity value of said range of velocity values with 
one or more corresponding ranges of time values. 



30 



72. 



The apparatus of claim 71 wherein said processor circuit is configured 
to, for each of a plurality of CMP gathers, define a window in which a 
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selected CMP location associated with said each CMP gather is 
centered. 

73. The apparatus of claim 72 wherein said processor circuit is configured 
5 to produce a semblance panel associated with said selected CMP 

location, in response to said range of velocity values and selected CMP 
gathers associated with CMP locations within said window. 

74. The apparatus of claim 73 wherein said processor circuit is configured 
10 to produce a time-velocity profile associated with said selected CMP 

location in response to said semblance panel and a set of velocity 
stacking weights. 

75. The apparatus of claim 74 wherein said processor circuit is configured 
15 to produce smoothing weights for each velocity value of said range of 

-velocity values and each time value of said one or more 

correspondence ranges of time values in response to respective 
products of total semblance and reciprocal velocity gradient associated 
with said each time value and said each velocity value. 

20 

76. The apparatus of claim 75 wherein said processor circuit is configured 
to produce a smooth time-velocity profile associated with said selected 
CMP location in response to said time-velocity profile and said 
smoothing weights. 

25 

77. The apparatus of claim 76 wherein said processor circuit is configured 
to produce laterally edited and laterally smoothed time-velocity profiles 
in response to a plurality of said smooth time-velocity profiles. 



30 



78. 



The apparatus of claim 70 wherein said processor circuit is configured 
to produce a migrated starting velocity field. 
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79. The apparatus of claim 78 wherein said processor circuit is configured 
to produce said migrated starting velocity field by, for said each CMP 
location, producing a normal moveout gather in response to said 
starting velocity field estimate and seismic data associated with said 

5 each CMP location. 

80. The apparatus of claim 79 wherein said processor circuit is configured 
to produce a CMP stack for said each CMP location in response to a 
corresponding said normal moveout gather. 

10 

81 . The apparatus of claim 80 wherein said processor circuit is configured 
to produce a simulated CMP gather for said each CMP location in 
response to a corresponding said CMP stack and a corresponding 
time-velocity profile associated with said each CMP location. 

15 

82. The apparatus of claim 81 wherein said processor circuit is configured 
to produce pre-stack imaged simulated gathers in response to a 
plurality of respective said simulated CMP gathers and respective said 
CMP locations. 

20 

83. The apparatus of claim 82 wherein said processor circuit is configured 
to produce normal moveout pre-stack imaged simulated gathers in 
response to said pre-stack imaged simulated gathers, said CMP 
locations and said starting velocity field estimate. 

25 

84. The apparatus of claim 83 wherein said processor circuit is configured 
to produce said migrated velocity by producing said migrated velocity 
field in response to said normal moveout pre-stack imaged simulated 
gathers and said starting velocity field estimate. 

30 

85. The apparatus of claim 84 wherein said processor circuit is configured 
to produce said migrated velocity field by, for said each CMP location 
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and for each of said common seismic events, finding a slope of a curve 
approximating a trend, relative to offset, in said time values, associated 
with said each common seismic event. 

The apparatus of claim 85 wherein said processor circuit is configured 
to find said slope of said curve by producing slant stack values for a 
plurality of time(tau)-slope pairs. 

The apparatus of claim 86 wherein said processor circuit is configured 
to, for said each CMP location, produce for each time value of said 
time(tau)-slope pairs a best slope value in response to said time(tau)- 
slope pairs and said slant stack values. 

The apparatus of claim 87 wherein said processor circuit is configured 
to produce a first temporary time-velocity representation in response to 
- said starting velocity field estimate and a plurality of said best slope 
values. 

The apparatus of claim 88 wherein said processor circuit is configured 
to produce said first temporary time-velocity representation by 
adjusting each velocity value at a given time t of a corresponding time- 
velocity representation of said starting velocity field estimate according 
to the relation: 
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where 



is an adjusted velocity value associated with said given 



time /; 



r 



P r 



is a Tau value equal to said given time t \ 
is said best slope value at r = t ; and 



Vi 



rep 



is said each velocity value at said given time / of said 



corresponding time-velocity representation. 



10 90. 
15 

91; 

20 

92. 

25 

93. 



The apparatus of claim 89 wherein said processor circuit is configured 
to condition said plurality of best slope values to produce conditioned 
best slope values and to produce an adjusted time-velocity 
representation in response to said temporary time-velocity 
representation and said conditioned best slope values according to 
saidTelation. " 

The apparatus of claim 90 wherein said processor circuit is configured 
to produce pre-stack imaged gathers by performing a 2-dimensional 
pre-stack imaging operation on seismic data associated with said 
control plane. 

The apparatus of claim 91 wherein said processor circuit is configured 
to produce normal moveout gathers in response to said migrated 
starting velocity field by performing respective normal moveout 
operations on said pre-stack imaged gathers. 

The apparatus of claim 92 wherein said processor circuit is configured 
to adjust said migrated starting velocity field in response to said normal 
moveout gathers and said migrated starting velocity field. 
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5 

95. 

10 

96. 



15 

97. 



20 

98. 



The apparatus of claim 93 wherein said processor circuit is configured 
to adjust said migrated starting velocity field by, for said each CMP 
location and for each of said common seismic events, finding a slope of 
a curve approximating a trend, relative to offset, in said time values, 
associated with said each common seismic event. 

The apparatus of claim 94 wherein said processor circuit is configured 
to find said slope of said curve by producing slant stack values for a 
plurality of time(tau)-slope pairs. 

The apparatus of claim 95 wherein said processor circuit is configured 
to, for said each CMP location, produce for each time value of said 
time(tau)-slope pairs a best slope value in response to said time(tau)- 
slope pairs and said slant stack values. 

The apparatus of claim 96 wherein said processor circuit is configured 
to produce a second temporary time-velocity representation in 
response to said migrated starting velocity field and a plurality of said 
best slope values. 

The apparatus of claim 97 wherein said processor circuit is configured 
to produce said second temporary time-velocity representation by 
adjusting each velocity value at a given time t of a corresponding time- 
velocity representation of said migrated starting velocity field according 
to the relation: 
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where 



is an adjusted velocity value associated with said given 



time /; 



P T 



T 



is a Tau value equal to said given time t \ 
is said best slope value at r = t\ and 



rep 



is said each velocity value at said given time / of said 



corresponding time-velocity representation. 



10 



15 



20 



99. The apparatus of claim 98 wherein said processor circuit is configured 
to condition said plurality of best slope values to produce conditioned 
best slope values and adjust said migrated starting velocity field in 
response to said temporary time-velocity representation and said 
conditioned best slope values according to said relation. 



100. The apparatus of claim 55 further comprising a processor circuit 
configured to determine whether the seismic data is 2-dimensional data 
or 3-dimensional data, determine whether an imaging velocity is 
desired to be produced, and perform a directional migration on the 
seismic data when the seismic data is 3-dimensional data and an 
imaging velocity is desired to be produced. 

101. The apparatus of claim 100 wherein said processor circuit is configured 
to determine whether or not a dip moveout velocity or said imaging 
velocity is desired to be produced, and if so, perform a dip moveout 
operation on the seismic data. 
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102. 



The apparatus of claim 101 wherein said processor circuit is configured 
to perform said dip moveout operation by performing a 2-dimensional 
dip moveout operation when the seismic data is said 2-dimensional 
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seismic data and perform a 3-dimensional dip moveout operation when 
the seismic data is 3-dimensional data. 

An apparatus for estimating a seismic velocity field from seismic data 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations 
aligned along edges of respective planes in a 3-dimensional region in 
the earth, the apparatus comprising: 

a) a component configured to designate one of said respective 
planes as a control plane; 

b) the apparatus of claim 55 for estimating a velocity field for said 
control plane; and 

c) a component configured to produce an optimized velocity field 
for one of said respective planes nearby said control plane by 
adjusting time-velocity values for each of the CMP locations 
associated with said control plane, in response to trends, 
relative to offset, in time values, associated with common 
seismic events, in seismic data associated with said one of said 
planes nearby said control plane until said optimized velocity 
field satisfies a condition. 

An apparatus for estimating a seismic velocity field from seismic data 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations 
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aligned along edges of respective planes in a 3-dimensional region in 
the earth, the apparatus comprising: 

a) a component configured to designate one of said respective 
planes as a control plane; 

b) the apparatus of claim 55 for estimating a velocity field for said 
control plane; and 

c) a component configured to produce an optimized velocity field 
for each of the respective planes in the 3-dimensional region by 
adjusting time-velocity values for each of the CMP locations 
associated with said control plane, in response to trends, 
relative to offset, in time values, associated with common 
seismic events, relative to offset, in the seismic data associated 
with said each of said respective planes until said optimized 
velocity field satisfies a condition. 

An apparatus for estimating a seismic velocity field from seismic data 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations 
aligned along edges of respective planes in a 3-dimensional region in 
the earth, the apparatus comprising: 

a) a component configured to designate one of the respective 
planes as a temporary control plane; 

b) the apparatus of claim 55 for estimating a velocity field for said 
temporary control plane; and 
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c) 

5 



10 

106. A method for producing an initial velocity field estimate for a control 
plane from seismic data associated with said control plane and 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
15 midpoint therebetween, the seismic data being arranged into common 

midpoint (CMP) gathers associated with respective CMP locations and 
the control plane having an edge intersecting a plurality of the CMP 
locations, the method comprising: 

20 a) producing a starting velocity field estimate from an initial range 

of velocity values and an initial range of time values; 

b) producing a migrated starting velocity field from said starting 
velocity field estimate and said seismic data; 

25 

c) producing pre-stack imaged gathers by performing a 2- 
dimensional pre-stack imaging process on said seismic data; 



a component configured to produce an optimized velocity field 
for successive ones of the respective planes in the 3- 
dimensional region by, for a given plane in the 3-dimensional 
region, adjusting time-velocity values of a predecessor one of 
the respective planes, starting with said temporary control 
plane, in response to trends, relative to offset, in time values, 
associated with common seismic events and a common CMP, 
in seismic data associated with said given plane until said 
optimized velocity field satisfies a condition. 



30 



producing normal moveout gathers in response to said migrated 
starting velocity field, including performing a normal moveout 
operation on said pre-stack imaged gathers; 
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adjusting said migrated starting velocity field in response to said 
normal moveout gathers and said migrated starting velocity field 
to produce a plurality of time-velocity values for each of the CMP 
locations, said plurality of said time-velocity values acting as 
said initial velocity field. 

107. The method of claim 106 wherein producing said starting velocity field 
estimate comprises defining a range of velocity values, defining one or 
more ranges of time values and associating each velocity value of said 
10 range of velocity values with one or more corresponding ranges of time 

values. 



e) 



5 



108. The method of claim 107 further comprising, for each of a plurality of 
CMP gathers, defining a window in which a selected CMP location 

15 associated with said each CMP gather is centered. 

109. The method of claim 108 further comprising producing a semblance 
panel associated with said selected CMP location, in response to said 
range of velocity values and selected CMP gathers associated with 

20 CMP locations within said window. 

110. The method of claim 109 further comprising producing a time-velocity 
profile associated with said selected CMP location in response to said 
semblance panel and a set of velocity stacking weights. 

25 

111. The method of claim 110 further comprising producing smoothing 
weights for each velocity value of said range of velocity values and 
each time value of said one or more corresponding ranges of time 
values in response to respective products of total semblance and 

30 reciprocal velocity gradient associated with said each time value and 

said each velocity value. 
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112. The method of claim 111 further comprising producing a smooth time- 
velocity profile associated with said selected CMP location in response 
to said time-velocity profile and said smoothing weights. 

5 113. The method of claim 112 further comprising producing laterally edited 

and laterally smoothed time-velocity profiles in response to a plurality 
of said smooth time-velocity profiles. 

114. The method of claim 106 wherein producing said migrated starting 
10 velocity field comprises, for said each of the CMP locations, producing 

a normal moveout gather in response to said starting velocity field 
estimate and seismic data associated with said each CMP locations. 

115. The method of claim 114 further comprising producing a CMP stack for 
15 said each CMP location in response to a corresponding said normal 

moveout gather. 

116. The method of claim 115 further comprising producing a simulated 
CMP gather for said each CMP location in response to a corresponding 

20 said CMP stack and a corresponding time-velocity profile associated 

with said each CMP location. 

117. The method of claim 116 further comprising producing pre-stack 
imaged simulated gathers in response to a plurality of respective said 

25 simulated CMP gathers and respective said CMP locations. 

118. The method of claim 117 further comprising producing normal moveout 
pre-stack imaged simulated gathers in response to said pre-stack 
imaged simulated gathers, said CMP locations and said starting 

30 velocity field estimate. 
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119. The method of claim 118 wherein producing said migrated velocity field 
comprises producing said migrated velocity field in response to said 
normal moveout pre-stack imaged simulated gathers and said starting 
velocity field estimate. 

5 

120. The method of claim 119 wherein producing said migrated velocity field 
comprises, for said each CMP location and for each of said common 
seismic events, finding a slope of a curve approximating a trend, 
relative to offset, in said time values, associated with said each 

1 0 common seismic event. 

121. The method of claim 120 wherein finding said slope of said curve 
comprises producing slant stack values for a plurality of time(tau)-slope 
pairs. 

15 

122. The method of claim 121 further comprising, for said each CMP 
location, producing for each time value of said time(tau)-slope pairs a 
best slope value in response to said time(tau)-slope pairs and said 
slant stack values. 

20 

123. The method of claim 122 further comprising producing a first temporary 
time-velocity representation in response to said starting velocity field 
estimate and a plurality of said best slope values. 

25 124. The method of claim 123 wherein producing said first temporary time- 

velocity representation comprises adjusting each velocity value at a 
given time t of a corresponding time-velocity representation of said 
starting velocity field estimate according to the relation: 
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Vadj(0 = 



10 




where 

VadjW is an adjusted velocity value associated with said given 
time t \ 

r is a Tau value equal to said given time / ; 

p T is said best slope value at r = t ; and 

V rep is saicl each velocit y value at said given time t of said 

corresponding time-velocity representation. 



125. The method of claim 124 further comprising conditioning said plurality 
of best slope values to produce conditioned best slope values and 
producing an adjusted time-velocity representation in response to said 
temporary time-velocity representation and said conditioned best slope 

15 values according to said relation. 

126. The method of claim 125 wherein adjusting said migrated starting 
velocity field comprises, for each of the CMP locations and for each of 
said common seismic events, finding a slope of a curve approximating 

20 said trend, relative to offset, in said time values, associated with said 

each common seismic event. 



25 



127. The method of claim 126 wherein finding said slope of said curve- 
comprises producing slant stack values for a plurality of time(tau)-slope 
pairs. 
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128. 



5 

129. 



10 130. 



15 • 



The method of claim 127 further comprising, for said each CMP 
location, producing for each time value of said time(tau)-slope pairs a 
best slope value in response to said time(tau)-slope pairs and said 
slant stack values. 

The method of claim 128 further comprising producing a second 
temporary time-velocity representation in response to said migrated 
starting velocity field and a plurality of said best slope values. 

The method of claim 129 wherein producing said second temporary 
time-velocity representation comprises adjusting each velocity value at 
a given time t of a corresponding time-velocity representation of said 
migrated starting velocity field according to a relation: 




where 



Vadj<» 



is an adjusted velocity value associated with said given 



P T 



T 



time t \ 

is a Tau value equal to said given time / ; 
is said best slope value at r = t ; and 



rep 



is said each velocity value at said given time t of said 



corresponding time-velocity representation. 



25 



131. 



The method of claim 130 further comprising conditioning said plurality 
of best slope values to produce conditioned best slope values and 
adjusting said migrated starting velocity field in response to said 
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temporary time-velocity representation and said conditioned best slope 
values according to said relation. 

132. A computer readable medium encoded with codes for directing a 
5 processor circuit to carry out the method of claim 106. 

133. A computer readable signal encoded with codes for directing a 
processor circuit to carry out the method of claim 106. 

10 134. An apparatus for producing an initial velocity field estimate for a control 

plane from seismic data associated with said control plane and 
comprising time-amplitude representations associated with source- 
receiver locations spaced apart by an offset distance and having a 
midpoint therebetween, the seismic data being arranged into common 

15 midpoint (CMP) gathers associated with respective CMP locations and 

the control plane having an edge- intersecting a plurality of the CMP 
locations, the apparatus comprising: 



a) means for producing a starting velocity field estimate from an 
20 initial range of velocity values and an initial range of time values; 

b) means for producing a migrated starting velocity field from said 
starting velocity field estimate and said seismic data; 

25 c) means for producing pre-stack imaged gathers by performing a 

2-dimensional pre-stack imaging process on said seismic data; 



30 



d) 



means for producing normal moveout gathers in response to 
said migrated starting velocity field, including performing a 
normal moveout operation on said pre-stack imaged gathers; 
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e) means for adjusting said migrated starting velocity field in 
response to said normal moveout gathers and said migrated 
starting velocity field to produce a plurality of time-velocity 
values for each of the CMP locations, said plurality of said time- 
5 velocity values acting as said initial velocity field. 

135. An apparatus for producing an initial velocity field estimate for a control 
plane from seismic data associated with said control plane and 
comprising time-amplitude representations associated with source- 
10 receiver locations spaced apart by an offset distance and having a 

midpoint therebetween, the seismic data being arranged into common 
midpoint (CMP) gathers associated with respective CMP locations and 
the control plane having an edge intersecting a plurality of the CMP 
locations, the apparatus comprising: 

15 

a) a component configured to produce a starting velocity field 
estimate from an initial range of velocity values and an initial 
range of time values; 

20 b) a component configured to produce a migrated starting velocity 

field from said starting velocity field estimate and said seismic 
data; 

c) a component configured to produce pre-stack imaged gathers 
25 by performing a 2-dimensional pre-stack imaging process on 

said seismic data; 

d) a component configured to produce normal moveout gathers in 
response to said migrated starting velocity field, including 

30 performing a normal moveout operation on said pre-stack 

imaged gathers; 
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e) a component configured to adjust said migrated starting velocity 
field in response to said normal moveout gathers and said 
migrated starting velocity field to produce a plurality of time- 
velocity values for each of the CMP locations, said plurality of 
5 said time-velocity values acting as said initial velocity field. 

136. The apparatus of claim 135 wherein said component configured to 
produce said starting velocity field estimate comprises a processor 
circuit configured to define a range of velocity values, define one or 

10 more ranges of time values and associate each velocity value of said 

range of velocity values with one or more corresponding ranges of time 
values. 

137. The apparatus of claim 136 wherein said processor circuit is configured 
15 to, for each of a plurality of CMP gathers, define a window in which a 

selected CMP -location associated with said each CMP gather is 
centered. 



138. The apparatus of claim 137 wherein said processor circuit is configured 
20 to produce a semblance panel associated with said selected CMP 

location, in response to said range of velocity values and selected CMP 
gathers associated with CMP locations within said window. 



139. The apparatus of claim 138 wherein said processor circuit is configured 
25 to produce a time-velocity profile associated with said selected CMP 

location in response to said semblance panel and a set of velocity 
stacking weights. 



140. The apparatus of claim 139 wherein said processor circuit is configured 
30 to produce smoothing weights for each velocity value of said range of 

velocity values and each time value of said one or more corresponding 
ranges of time values in response to respective products of total 
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semblance and reciprocal velocity gradient associated with said each 
time value and said each velocity value. 

141. The apparatus of claim 140 wherein said processor circuit is configured 
5 to produce a smooth time-velocity profile associated with said selected 

CMP location in response to said time-velocity profile and said 
smoothing weights. 

142. The apparatus of claim 141 wherein said processor circuit is configured 
10 to produce laterally edited and laterally smoothed time-velocity profiles 

in response to a plurality of said smooth time-velocity profiles. 

143. The apparatus of claim 135 wherein said component configured to 
produce said migrated starting velocity field comprises said processor 

15 and wherein said processor is configured to, for said each of the CMP 

location, produce a normal moveout gather in response to said starting 
velocity field estimate and seismic data associated with said each CMP 
location. 

20 144. The apparatus of claim 143 wherein said processor circuit is configured 

to produce a CMP stack for said each CMP location in response to a 
corresponding said normal moveout gather. 

145. The apparatus of claim 144 wherein said processor circuit is configured 
25 to produce a simulated CMP gather for said each CMP location in 

response to a corresponding said CMP stack and a corresponding 
time-velocity profile associated with said each CMP location. 

146. The apparatus of claim 145 wherein said processor circuit is configured 
30 to produce pre-stack imaged simulated gathers in response to a 

plurality of respective said simulated CMP gathers and respective said 
CMP locations. 
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147. The apparatus of claim 146 wherein said processor circuit is configured 
to produce normal moveout pre-stack imaged simulated gathers in 
response to said pre-stack imaged simulated gathers, said CMP 
locations and said starting velocity field estimate. 

5 

148. The apparatus of claim 147 wherein said processor circuit is configured 
to produce said migrated velocity field by producing said migrated 
velocity field in response to said normal moveout pre-stack imaged 
simulated gathers and said starting velocity field estimate. 

10 

149. The apparatus of claim 148 wherein said processor circuit is configured 
to produce said migrated velocity field by, for said each CMP location 
and for each of said common seismic events, finding a slope of a curve 
approximating a trend, relative to offset, in said time values, associated 

1 5 with said each common seismic event. 

150. The apparatus of claim 149 wherein said processor circuit is configured 
to find said slope of said curve by producing slant stack values for a 
plurality of time(tau)-slope pairs. 

20 

151. The apparatus of claim 150 wherein said processor circuit is configured 
to, for said each CMP location, produce for each time value of said 
time(tau)-slope pairs a best slope value in response to said time(tau)- 
slope pairs and said slant stack values. 

25 

152. The apparatus of claim 151 wherein said processor circuit is configured 
to produce a first temporary time-velocity representation in response to 
said starting velocity field estimate and a plurality of said best slope 
values. 

30 

153. The apparatus of claim 152 wherein said processor circuit is configured 
to produce said first temporary time-velocity representation by 
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adjusting each velocity value at a given time t of a corresponding time- 
velocity representation of said starting velocity field estimate according 
to the relation: 



VadjW 




where 



V a dj(0 is an adjusted velocity value associated with said given 

10 time t\ 

r is a Tau value equal to said given time / ; 

p T is said best slope value at r = t ; and 

V rep ' s sa ' c ' eac ^ velocity value at said given time t of said 

corresponding time-velocity representation. 

15 

154. The apparatus of claim 153 wherein said processor circuit is configured 
to condition said plurality of best slope values to produce conditioned 
best slope values and produce an adjusted time-velocity representation 
in response to said temporary time-velocity representation and said 
20 conditioned best slope values according to said relation. 



155. The apparatus of claim 154 wherein said component configured to 
adjust said migrated starting velocity field comprises said processor 
circuit configured to, for each of the CMP locations and for each of said 
25 common seismic events, find a slope of a curve approximating said 

trend, relative to offset, in said time values, associated with said each 
common seismic event. 
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10 



15 



20 



156. The apparatus of claim 155 wherein said processor circuit is configured 
to find said slope of said curve by producing slant stack values for a 
plurality of time(tau)-slope pairs. 

157. The apparatus of claim 156 wherein said processor circuit is configured 
to, for said each CMP location, produce for each time value of said 
time(tau)-slope pairs a best slope value in response to said time(tau)- 
slope pairs and said slant stack values. 

158. The apparatus of claim 157 wherein said processor circuit is configured 
to produce a second temporary time-velocity representation in 
response to said migrated starting velocity field and a plurality of said 
best slope values. 

159. The apparatus of claim 158 wherein said processor circuit is configured 
to produce said second temporary time-velocity representation by 
adjusting each velocity value at a given time t of a corresponding time- 
velocity representation of said migrated starting velocity field according 
to a relation: 




where 



is an adjusted velocity value associated with said given 



T 



P r 



time /; 

is a Tau value equal to said given time t\ 
is said best slope value at r = t\ and 
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' rep is saicl each ve,oc 'ty value at said given time / of said 

corresponding time-velocity representation. 



160. The apparatus of claim 159 wherein said processor circuit is configured 
5 to condition said plurality of best slope values to produce conditioned 

best slope values and adjust said migrated starting velocity field in 
response to said temporary time-velocity representation and said 
conditioned best slope values according to said relation. 

10 161. A method for producing a starting velocity field estimate from seismic 

data comprising time-amplitude representations associated with 
source-receiver locations spaced apart by an offset distance and 
having a midpoint therebetween, the seismic data being arranged into 
common midpoint (CMP) gathers associated with respective CMP 

15 locations, the method comprising: 

a) defining a range of velocity values, defining one or more ranges 
of time values and associating each velocity value of said range 
of velocity values with one or more corresponding ranges of time 

20 values; and 

b) for each of a plurality of CMP gathers: 

i) defining a window in which a selected CMP location 
25 associated with said each CMP gather is centered; 

ii) producing a semblance panel associated with said 
selected CMP location, in response to said range of 
velocity values and selected CMP gathers associated 

30 with CMP locations within said window; and 
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iii) producing a time-velocity profile associated with said 
selected CMP location in response to said semblance 
panel and a set of velocity stacking weights; 

5 iv) producing smoothing weights for each velocity value of 

said range of velocity values and each time value of said 
one or more ranges of time values in response to 
respective products of total semblance and reciprocal 
velocity gradient associated with said each time value 
10 and said each velocity value; and 

v) producing a smooth time-velocity profile associated with 
said selected CMP location in response to said time- 
velocity profile and said smoothing weights. 

15 

162. The method of claim 161 further comprising producing laterally edited 
and laterally smoothed time-velocity profiles in response to a plurality 
of said smooth time-velocity profiles. 

20 163. The method of claim 161 further comprising selecting CMP gathers of 

the seismic data that meet selection criteria based on CMP fold and 
offset distribution, for use in defining said window. 

164. A computer readable medium encoded with codes for directing a 
25 processor circuit to carry out the method of claim 161 . 

165. A computer readable signal encoded with codes for directing a 
processor circuit to carry out the method of claim 161. 



30 166. 



An apparatus for producing a starting velocity field estimate from 
seismic data comprising time-amplitude representations associated 
with source-receiver locations spaced apart by an offset distance and 
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having a midpoint therebetween, the seismic data being arranged into 
common midpoint (CMP) gathers associated with respective CMP 
locations, the apparatus comprising: 



a) means for defining a range of velocity values, defining one or 
more ranges of time values and associating each velocity value 
of said range of velocity values with one or more corresponding 
ranges of time values; and 

b) means for, for each of a plurality of CMP gathers: 

i) defining a window in which a selected CMP location 
associated with said each CMP gather is centered; 

ii) producing a semblance panel associated with said 
selected CMP location, in response to said range of 
velocity values and selected CMP gathers associated 
with CMP locations within said window; and 

iii) producing a time-velocity profile associated with said 
selected CMP location in response to said semblance 
panel and a set of velocity stacking weights; 

iv) producing smoothing weights for each velocity value of 
said range of velocity values and each time value of said 
one or more ranges of time values in response to 
respective products of total semblance and reciprocal 
velocity gradient associated with said each time value 
and said each velocity value; and 



-126- 

v) producing a smooth time-velocity profile associated with 
said selected CMP location in response to said time- 
velocity profile and said smoothing weights. 

5 167. An apparatus for producing a starting velocity field estimate from 

seismic data comprising time-amplitude representations associated 
with source-receiver locations spaced apart by an offset distance and 
having a midpoint therebetween, the seismic data being arranged into 
common midpoint (CMP) gathers associated with respective CMP 
10 locations, the apparatus comprising: 

a) a component configured to define a range of velocity values, 
define one or more ranges of time values and associate each 
velocity value of said range of velocity values with one or more 

1 5 corresponding ranges of time values; and 

b) a component configured to for each of a plurality of CMP 
gathers: 

20 i) define a window in which a selected CMP location 

associated with said each CMP gather is centered; 

ii) produce a semblance panel associated with said selected 
CMP location, in response to said range of velocity 
25 values and selected CMP gathers associated with CMP 

locations within said window; and 



30 



iii) 



produce a time-velocity profile associated with said 
selected CMP location in response to said semblance 
panel and a set of velocity stacking weights; 
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iv) produce smoothing weights for each velocity value of said 
range of velocity values and each time value of said one 
or more ranges of time values in response to respective 
products of total semblance and reciprocal velocity 
5 gradient associated with said each time value and said 

each velocity value; and 



v) produce a smooth time-velocity profile associated with 
said selected CMP location in response to said time- 
10 velocity profile and said smoothing weights. 

168. The apparatus of claim 167 further comprising a component configured 
to produce laterally edited and laterally smoothed time-velocity profiles 
in response to a plurality of said smooth time-velocity profiles. 

15 

169. The apparatus of claim 167 wherein said component configured to 
define a window is configured to select CMP gathers of the seismic 
data that meet selection criteria based on CMP fold and offset 
distribution, for use in defining said window. 

20 

170. A method for producing a migrated velocity field in response to seismic 
data comprising time-amplitude representations associated with 
source-receiver locations spaced apart by an offset distance and 
having a midpoint therebetween, the seismic data being arranged into 

25 common midpoint (CMP) gathers associated with respective CMP 

locations and a starting velocity field estimate comprising time-velocity 
profiles associated with respective said CMP locations, the method 
comprising: 



30 



for each of the CMP locations: 
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i) producing a normal moveout gather in response to the 
starting velocity field estimate and the seismic data 
associated with said each CMP location; 

5 ii) producing a CMP stack in response to a corresponding 

said normal moveout gather; and 

Hi) producing a simulated CMP gather in response to a 
corresponding said CMP stack and a corresponding time- 
10 velocity profile; 

b) producing pre-stack imaged simulated gathers in response to a 
plurality of respective said simulated CMP gathers and 
respective said CMP locations; 

15 

c) producing normal moveout pre-stack imaged simulated gathers 
in response to said pre-stack imaged simulated gathers, said 
CMP locations and the starting velocity field estimate; and 

20 d) producing the migrated velocity field in response to said normal 

moveout pre-stack imaged simulated gathers and the starting 
velocity field estimate. 

171. The method of claim 170 wherein producing the migrated velocity field 
25 comprises, for said each CMP location and for each of a plurality of 

common seismic events, finding a slope of a curve approximating a 
trend, relative to offset, in time values, associated with said each 
common seismic event. 



30 



172. The method of claim 171 wherein finding said slope of said curve 
comprises producing slant stack values for a plurality of time(tau)-slope 
pairs. 
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173. The method of claim 172 further comprising, for said each CMP 
location, producing for each time value of said time(tau)-slope pairs a 
best slope value in response to said time(tau)-slope pairs and said 
slant stack values. 

5 

174. The method of claim 173 further comprising producing a temporary 
time-velocity representation in response to the starting velocity field 
estimate and a plurality of said best slope values. 

10 175. The method of claim 174 wherein producing said temporary time- 

velocity representation comprises adjusting each velocity value at a 
given time t of a corresponding time-velocity representation of the 
starting velocity field estimate according to a relation: 



15 



Vadj«) = 



V 



rep 



where 



V ac lj(t) is an adjusted velocity value associated with said given 

20 time t\ 

r is a Tau value equal to said given time / ; 

p T is said best slope value at r = / ; and 

V rep is said each ve,ocit y value at said given time t of said 

corresponding time-velocity representation. 



25 
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176. The method of claim 175 further comprising conditioning said plurality 
of best slope values to produce conditioned best slope values and 
producing an adjusted time-velocity representation in response to said 
temporary time-velocity representation and said conditioned best slope 

5 values according to said relation. 

177. A computer readable medium encoded with codes for directing a 
processor circuit to carry out the method of claim 170. 

10 178. A computer readable signal encoded with codes for directing a 

processor circuit to carry out the method of claim 170. 



179. An apparatus for producing a migrated velocity field in response to 
seismic data comprising time-amplitude representations associated 

15 with source-receiver locations spaced apart by an offset distance and 

having a midpoint therebetween, the seismic data being arranged into 
common midpoint (CMP) gathers associated with respective CMP 
locations and a starting velocity field estimate comprising time-velocity 
profiles associated with respective said CMP locations, the apparatus 

20 comprising: 

a) means for, for each of the CMP locations: 



i) producing a normal moveout gather in response to the 
25 starting velocity field estimate and the seismic data 

associated with said each CMP location; 

ii) producing a CMP stack in response to a corresponding 
said normal moveout gather; and 

30 
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iii) producing a simulated CMP gather in response to a 
corresponding said CMP stack and a corresponding time- 
velocity profile; 

5 b) means for producing pre-stack imaged simulated gathers in 

response to a plurality of respective said simulated CMP gathers 
and respective said CMP locations; 

c) means for producing normal moveout pre-stack imaged 
10 simulated gathers in response to said pre-stack imaged 

simulated gathers, said CMP locations and the starting velocity 
field estimate; and 

d) means for producing the migrated velocity field in response to 
15 said normal moveout pre-stack imaged simulated gathers and 

the starting velocity field estimate. 

180. An apparatus for producing a migrated velocity field in response to 
seismic data comprising time-amplitude representations associated 

20 with source-receiver locations spaced apart by an offset distance and 

having a midpoint therebetween, the seismic data being arranged into 
common midpoint (CMP) gathers associated with respective CMP 
locations and a starting velocity field estimate comprising time-velocity 
profiles associated with respective said CMP locations, the apparatus 

25 comprising: 

a) a component configured to, for each of the CMP locations: 



30 



produce a normal moveout gather in response to the 
starting velocity field estimate and the seismic data 
associated with said each CMP location; 
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ii) produce a CMP stack in response to a corresponding 
said normal moveout gather; and 

iii) produce a simulated CMP gather in response to a 
5 corresponding said CMP stack and a corresponding time- 
velocity profile; 

b) a component configured to produce pre-stack imaged simulated 
gathers in response to a plurality of respective said simulated 

10 CMP gathers and respective said CMP locations; 

c) a component configured to produce normal moveout pre-stack 
imaged simulated gathers in response to said pre-stack imaged 
simulated gathers, said CMP locations and the starting velocity 

15 field estimate; and 

d) a component configured to produce the migrated velocity field in 
response to said normal moveout pre-stack imaged simulated 
gathers and the starting velocity field estimate. 

20 

181. The apparatus of claim 180 wherein said component configured to 
produce the migrated velocity field comprises a processor circuit 
configured to, for said each CMP location and for each of a plurality of 
common seismic events, find a slope of a curve approximating a trend, 
25 relative to offset, in time values, associated with said each common 

seismic event. 



30 



182. The apparatus of claim 181 wherein said processor circuit is configured 
to find said slope of said curve by producing slant stack values for a 
plurality of time(tau)-slope pairs. 



10 



-133- 

183. The apparatus of claim 182 wherein said processor circuit is configured 
to, for said each CMP location, produce for each time value of said 
time(tau)-slope pairs a best slope value in response to said time(tau)- 
slope pairs and said slant stack values. 

184. The apparatus of claim 183 wherein said processor circuit is configured 
to produce a temporary time-velocity representation in response to the 
starting velocity field estimate and a plurality of said best slope values. 

185. The apparatus of claim 184 wherein said processor circuit is configured 
to produce said temporary time-velocity representation by adjusting 
each velocity value at a given time t of a corresponding time-velocity 
representation of the starting velocity field estimate according to a 
relation: 



15 



Vadj(0 = 




20 



25 



where 

T 

P r 
V rep 



is an adjusted velocity value associated with said given 
time t\ 

is a Tau value equal to said given time t\ 
is said best slope value at r ~t \ and 

is said each velocity value at said given time t of said 

corresponding time-velocity representation. 
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186. The apparatus of claim 185 where said processor circuit is configured 
to condition said plurality of best slope values to produce conditioned 
best slope values and produce an adjusted time-velocity representation 
in response to said temporary time-velocity representation and said 
5 conditioned best slope values according to said relation. 



187. A method for producing an output velocity field in response to an input 
velocity field and seismic data comprising time-amplitude 
representations associated with source-receiver locations spaced apart 
10 by an offset distance and having a midpoint therebetween, the seismic 

data being arranged into common midpoint (CMP) gathers associated 
with respective CMP locations, the method comprising: 



a) for each of the CMP locations and for each of a plurality of 
15 common seismic events: 

i) finding a slope of a curve approximating a trend, relative 
to offset, in time values, associated with said each 
common seismic event, by producing slant stack values 

20 for a plurality of time(tau)-slope pairs; 

ii) producing for each time value of said time(tau)-slope 
pairs a best slope value in response to said time(tau)- 
slope pairs and said slant stack values; 



25 



iii) producing a temporary time-velocity representation in 
response to the input velocity field and a plurality of said 
best slope values; and 



30 



iv) 



conditioning said plurality of best slope values to produce 
conditioned best slope values and producing an adjusted 
time-velocity representation in response to said 
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temporary time-velocity representation and said 
conditioned best slope values. 



10 



188. The method of claim 187 wherein producing said temporary time- 
velocity representation comprises adjusting each velocity value at a 
given time t of a corresponding time-velocity representation of the input 
velocity field according to a relation: 



VadjM= ! 



i 



Vrep 



where 



VadjW is an ac| j usted velocity value associated with said given 
time t \ 

15 r is a Tau value equal to said given time t ; 

p T is said best slope value at r = t ; and 

V rep is said each velocity value at said given time t of said 

corresponding time-velocity representation. 

20 189. The method of claim 187 wherein producing said slant stack values 

comprises storing said slant stack values in a buffer as said slant stack 
values are produced. 



190. 

25 



The method of claim 189 wherein producing said slant stack values 
comprises producing said slant stack values until said buffer reaches a 
fill level. 
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191. The method of claim 190 further comprising producing edited and 
smoothed tau-slope values for said each CMP location. 

192. The method of claim 191 wherein producing edited and smoothed tau- 
slope values comprises performing a weighted lateral median filtering 
on said time(tau)-slope pairs. 

193. The method of claim 187 further comprising producing a weight value 
for each velocity value of said temporary time-velocity representation. 

194. The method of claim 193 further comprising producing a best slope 
reliability value for each of said plurality of best slope values. 



195. The method of claim 194 further comprising producing said conditioned 
15 best slope values in response to said best slope values, said best 

slope reliability values and said weight values. 

196. The method of claim 188 wherein producing said adjusted time-velocity 
representation comprises producing said adjusted time-velocity 

20 representation in response to said temporary time-velocity 

representation and said conditioned best slope values according to 
said relation, using said temporary time velocity representation as said 
corresponding time-velocity representation. 

25 197. A computer readable medium encoded with codes for directing a 

processor circuit to carry out the method of claim 187. 

198. A computer readable signal encoded with codes for directing a 
processor circuit to carry out the method of claim 187. 



30 



199. An apparatus for producing an output velocity field in response to an 
input velocity field and seismic data comprising time-amplitude 
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representations associated with source-receiver locations spaced apart 
by an offset distance and having a midpoint therebetween, the seismic 
data being arranged into common midpoint (CMP) gathers associated 
with respective CMP locations, the apparatus comprising: 

5 

a) for each of the CMP locations and for each of a plurality of 
common seismic events: 

i) means for finding a slope of a curve approximating a 
10 trend, relative to offset, in time values, associated with 

said each common seismic event, by producing slant 
stack values for a plurality of time(tau)-slope pairs; 

ii) means for producing for each time value of said 
15 time(tau)-slope pairs a best slope value in response to 

said time(tau)-slope pairs and said slant stack values; 

iii) means for producing a temporary time-velocity 
representation in response to the input velocity field and a 

20 plurality of said best slope values; and 

iv) means for conditioning said plurality of best slope values 
to produce conditioned best slope values and producing 
an adjusted time-velocity representation in response to 

25 said temporary time-velocity representation and said 

conditioned best slope values, for each of the CMP 
locations and for each of a plurality of common seismic 
events. 



30 



200. 



An apparatus for producing an output velocity field in response to an 
input velocity field and seismic data comprising time-amplitude 
representations associated with source-receiver locations spaced apart 
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by an offset distance and having a midpoint therebetween, the seismic 
data being arranged into common midpoint (CMP) gathers associated 
with respective CMP locations, the apparatus comprising: 



5 a) a processor circuit configured to for each of the CMP locations 

and for each of a plurality of common seismic events: 

i) find a slope of a curve approximating a trend, relative to 
offset, in time values, associated with said each common 

10 seismic event, by producing slant stack values for a 

plurality of time(tau)-slope pairs; 

ii) produce for each time value of said time(tau)-slope pairs 
a best slope value in response to said time(tau)-slope 

1 5 pairs and said slant stack values; 

iii) produce a temporary time-velocity representation in 
response to the input velocity field and a plurality of said 
best slope values; and 

20 

iv) condition said plurality of best slope values to produce 
conditioned best slope values and produce an adjusted 
time-velocity representation in response to said 
temporary time-velocity representation and said 

25 conditioned best slope values. 

201. The apparatus of claim 200 wherein said processor circuit is configured 
to produce said temporary time-velocity representation by adjusting 
each velocity value at a given time t of a corresponding time-velocity 
30 representation of the input velocity field according to a relation: 
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10 




where 

VadjW is an ac, j ustecl velocity value associated with said given 
time t\ 

r is a Tau value equal to said given time t ; 

p T is said best slope value at r = t ; and 

V rep is said eacl1 velocity value at said given time t of said 

corresponding time-velocity representation. 

202. The apparatus of claim 200 wherein said processor circuit is configured 
to produce said slant stack values by storing said slant stack values in 
a buffer as said slant stack values are produced. 



15 203. The apparatus of claim 202 wherein said processor circuit is configured 

to produce said slant stack values by producing said slant stack values 
until said buffer reaches a fill level. 

204. The apparatus of claim 203 wherein said processor circuit is configured 
20 to produce edited and smoothed tau-slope values for said each CMP 

location. 



25 



205. 



The apparatus of claim 204 wherein said processor circuit is configured 
to produce edited and smoothed tau-slope values by performing a 
weighted lateral median filtering on said time(tau)-slope pairs. 
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206. The apparatus of claim 200 wherein said processor circuit is configured 
to produce a weight value for each velocity value of said temporary 
time-velocity representation. 

5 207. The apparatus of claim 206 wherein said processor circuit is configured 

to produce a best slope reliability value for each of said plurality of best 
slope values. 

208. The apparatus of claim 207 wherein said processor circuit is configured 
10 to produce said conditioned best slope values in response to said best 

slope values, said best slope reliability values and said weight values. 

209. The apparatus of claim 201 wherein said processor circuit is configured 
to produce said adjusted time-velocity representation by producing said 

15 adjusted time-velocity representation in response to said temporary 

time-velocity representation and said conditioned best slope values 
according to said relation, and use said temporary time velocity 
representation as said corresponding time-velocity representation. 

20 210. A method for producing an optimized velocity field in response to an 

initial velocity field and seismic data comprising time-amplitude 
representations associated with source-receiver locations spaced apart 
by an offset distance and having a midpoint therebetween, the seismic 
data being arranged into common midpoint (CMP) gathers associated 

25 with respective CMP locations, the method comprising: 

a) defining a temporary initial velocity field associated with said 
seismic data; 

30 b) producing a set of gathers in response to said temporary initial 

velocity field, according to a type of velocity field desired to be 
produced. 
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c) producing an adjusted velocity field in response to said set of 
gathers and said temporary initial velocity field. 

5 d) outputting said adjusted velocity field when said adjusted 

velocity field meets optimization criteria. 

211. The method of claim 210 further comprising producing a new adjusted 
velocity field using said adjusted velocity field as a new temporary 

10 initial velocity field when said optimization criteria are not met. 

212. The method of claim 210 further comprising identifying said adjusted 
velocity field as said optimized velocity field when said optimization 
criteria are met. 

15 

213. The method of claim 210 wherein producing said set of gathers 
comprises producing sets of common offset migration image gathers in 
response to the seismic data and said initial velocity field when an 
imaging velocity field is desired to be produced. 

20 

214. The method of claim 213 wherein producing said set of gathers 
comprises performing a normal moveout operation on the seismic data 
when said imaging velocity field is not desired to be produced. 

25 215. The method of claim 210 wherein outputting said adjusted velocity field 

when said adjusted velocity field meets optimization criteria comprises 
determining whether velocity values of said adjusted velocity field are 
within a range of corresponding velocity values of said temporary initial 
velocity field. 

30 

216. The method of claim 210 wherein producing said adjusted velocity field 
comprises, for each of the CMP locations and for each of a plurality of 
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5 217. 



218. 

10 



219. 

15 

220. 

20 



common seismic events, finding a slope of a curve approximating a 
trend, relative to offset, in time values, associated with said each 
common seismic event. 

The method of claim 216 wherein finding said slope of said curve 
comprises producing slant stack values for a plurality of time(tau)-slope 
pairs. 

The method of claim 217 further comprising, for said each CMP 
location, producing for each time value of said time(tau)-slope pairs a 
best slope value in response to said time(tau)-slope pairs and said 
slant stack values. 

The method of claim 218 further comprising producing a temporary 
time- velocity representation in response to the initial velocity field and a 
plurality of said best slope values. 

The method of claim 219 wherein producing said temporary time- 
velocity representation comprises adjusting each velocity value at a 
given time t of a corresponding time-velocity representation of the initial 
velocity field according to a relation: 




where 



V a( jj(t) is an adjusted velocity value associated with said given 
time /; 
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P T 



T 



is a Tau value equal to said given time /; 
is said best slope value at r = / ; and 



rep 



is said each velocity value at said given time t of said 



corresponding time-velocity representation. 



5 

221. 

10 

222. 

15 

223. 

20 224. 

225. 

25 

226. 

30 227. 



The method of claim 220 further comprising conditioning said plurality 
of best slope values to produce conditioned best slope values and 
producing an adjusted time-velocity representation in response to said 
temporary time-velocity representation and said conditioned best slope 
values according to said relation. 

The method of claim 217 wherein producing slant stack values 
comprises storing said slant stack values in a buffer as said slant stack 
values are produced. 

The method of claim 222 wherein producing said slant stack values 
comprises producing said slant stack values until said buffer reaches a 
fill level. 

The method of claim 223 further comprising producing edited and 
smoothed tau-slope values for said each CMP location. 

The method of claim 224 wherein producing edited and smoothed tau- 
slope values comprises performing a weighted lateral median filtering 
on said time(tau)-slope pairs. 

The method of claim 225 further comprising producing a weight value 
for each velocity value of said temporary time-velocity representation. 

The method of claim 226 further comprising producing a best slope 
reliability value for each of said plurality of best slope values. 
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228. The method of claim 227 further comprising producing said conditioned 
best slope values in response to said plurality of best slope values, a 
plurality of said best slope reliability values and a plurality of said 
weight values. 

5 

229. A computer readable medium encoded with codes for directing a 
processor circuit to carry out the method of claim 210. 

230. A computer readable signal encoded with codes for directing a 
10 processor circuit to carry out the method of claim 210. 

231. An apparatus for producing an optimized velocity field in response to 
an initial velocity field and seismic data comprising time-amplitude 
representations associated with source-receiver locations spaced apart 

15 by an offset distance and having a midpoint therebetween, the seismic 

data being arranged into common midpoint (CMP) gathers associated 
with respective CMP locations, the apparatus comprising: 

a) means for defining a temporary initial velocity field associated 
20 with said seismic data; 

b) means for producing a set of gathers in response to said 
temporary initial velocity field, according to a type of velocity 
field desired to be produced. 



25 



c) means for producing an adjusted velocity field in response to 
said set of gathers and said temporary initial velocity field. 



30 



d) 



means for outputting said adjusted velocity field when said 
adjusted velocity field meets optimization criteria. 
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10 



232. An apparatus for producing an optimized velocity field in response to 
an initial velocity field and seismic data comprising time-amplitude 
representations associated with source-receiver locations spaced apart 
by an offset distance and having a midpoint therebetween, the seismic 
data being arranged into common midpoint (CMP) gathers associated 
with respective CMP locations, the apparatus comprising: 

a) a component configured to define a temporary initial velocity 
field associated with said seismic data; 

b) a component configured to produce a set of gathers in response 
to said temporary initial velocity field, according to a type of 
velocity field desired to be produced; 



15 c) a component configured to produce an adjusted velocity field in 

response to said set of gathers and said temporary initial 
velocity field. 

d) a component configured to output said adjusted velocity field 
20 when said adjusted velocity field meets optimization criteria. 



233. The apparatus of claim 232 wherein said component configured to 
produce an adjusted velocity field is configured to produce a new 
adjusted velocity field using said adjusted velocity field as a new 

25 temporary initial velocity field when said optimization criteria are not 

met. 

234. The apparatus of claim 232 wherein said component configured to 
output said adjusted velocity field is configured to identify said adjusted 

30 velocity field as said optimized velocity field when said optimization 

criteria are met. 
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235. The apparatus of claim 232 wherein said component configured to 
produce said set of gathers is configured to produce sets of common 
offset migration image gathers in response to the seismic data and said 
initial velocity field when an imaging velocity field is desired to be 

5 produced. 

236. The apparatus of claim 235 wherein said component configured to 
produce said set of gathers is configured to perform a normal moveout 
operation on the seismic data when said imaging velocity field is not 

10 desired to be produced. 

237. The apparatus of claim 232 wherein said component configured to 
output said adjusted velocity field when said adjusted velocity field 
meets optimization criteria is configured to determine whether velocity 

15 values of said adjusted velocity field are within a range of 

corresponding velocity values of said temporary initial velocity field. 

238. The apparatus of claim 232 wherein said component configured to 
produce said adjusted velocity field includes a processor circuit 

20 configured to, for each of the CMP locations and for each of a plurality 

of common seismic events, find a slope of a curve approximating a 
trend, relative to offset, in time values, associated with said each 
common seismic event. 

25 239. The apparatus of claim 238 wherein said processor circuit is configured 

to find said slope of said curve by producing slant stack values for a 
plurality of time(tau)-slope pairs. 



240. The apparatus of claim 239 wherein said processor circuit is configured 
30 to, for said each CMP location, produce for each time value of said 

time(tau)-slope pairs a best slope value in response to said time(tau)- 
slope pairs and said slant stack values. 
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241 . The apparatus of claim 240 wherein said processor circuit is configured 
to produce a temporary time-velocity representation in response to the 
initial velocity field and a plurality of said best slope values. 

242. The apparatus of claim 241 wherein said processor circuit is configured 
to produce said temporary time-velocity representation by adjusting 
each velocity value at a given time t of a corresponding time-velocity 
representation of the initial velocity field according to a relation: 



VadjW 



1 



where 



VadjW is an ad j usted velocity value associated with said given 
time t\ 

15 r is a Tau value equal to said given time t ; 

p r is said best slope value at r = t ; and 

V rep is saicl each velocity value at said given time t of said 

corresponding time-velocity representation. 



20 243. The apparatus of claim 242 wherein said processor circuit is configured 

to condition said plurality of best slope values to produce conditioned 
best slope values and to produce an adjusted time-velocity 
representation in response to said temporary time-velocity 
representation and said conditioned best slope values according to 

25 said relation. 
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244. The apparatus of claim 239 wherein said processor circuit is configured 
to produce slant stack values by storing said slant stack values in a 
buffer as said slant stack values are produced. 

5 245. The apparatus of claim 244 wherein said processor circuit is configured 

to produce said slant stack values by producing said slant stack values 
until said buffer reaches a fill level. 

246. The apparatus of claim 245 wherein said processor circuit is configured 
10 to produce edited and smoothed tau-slope values for said each CMP 

location. 

247. The apparatus of claim 246 wherein said processor circuit is configured 
to produce edited and smoothed tau-slope values by performing a 

15 weighted lateral median filtering on said time(tau)-slope pairs. 

248. The apparatus of claim 247 wherein said processor circuit is configured 
to produce a weight value for each velocity value of said temporary 
time-velocity representation. 

20 

249. The apparatus of claim 248 wherein said processor circuit is configured 
to produce a best slope reliability value for each of said plurality of best 
slope values. 

25 250. The apparatus of claim 249 wherein said processor circuit is configured 

to produce said conditioned best slope values in response to said 
plurality of best slope values, a plurality of said best slope reliability 
values and a plurality of said weight values. 



